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O ■ Abstract 

(N 

The discrete variational identity under general bilinear forms on semi-direct sums 
of Lie algebras is established. The constant 7 involved in the variational identity is 
determined through the corresponding solution to the stationary discrete zero curvature 
equation. An application of the resulting variational identity to a class of semi-direct 
sums of Lie algebras in the Volterra lattice case furnishes Hamiltonian structures for 
the associated integrable couplings of the Volterra lattice hierarchy. 
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1 Introduction 

An algebraic approach to integrable couplings [H E] was recently presented, based on the 
concept of semi-direct sums of Lie algebras j3l 0]. There exist plenty of examples of both 
continuous and discrete integrable couplings belonging to such a class of integrable equa- 
tions [1]-[T0]. The corresponding results show various mathematical structures that integrable 
equations possess, and provide a powerful tool to analyze integrable equations, particularly, 
multi-component integrable equations and integrable couplings [HI E] . Observe that a gen- 
eral Lie algebra can be decomposed into a semi-direct sum of a solvable Lie algebra and 
a semisimple Lie algebra [12]. The semi-direct sum decomposition of Lie algebras allows 
for more classifications of integrable equations supplementing existing theories [T3im] . e.g 
classifications within the areas of symmetry reductions P2J [TU] and Lax pairs [T7] . 

Let G be a matrix loop algebra, E be the shift operator, and D denote the forward 
difference operator E — 1, i.e., D = E — 1. Traditionally, we write 

(E m f)(n) = f {m \n) = f(n + m), m,ne Z, (1.1) 
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and define an inverse of the difference operator E — E l as follows (see, e.g., [T£] ' 
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((E-E)- 1 f)(n) = -{ Y, f(n+l + 2m)-J2f(n-l + 2m)), n e Z, (1.2) 

— oo rra=l 

where / is an expression depending on the lattice variable n. The corresponding inverses of 
the forward and backward difference operators are determined by 

(E - I)' 1 = (E- E- 1 )- 1 ^ + E- 1 ), (1 - E- 1 )- 1 = (E- E- 1 )-\E + 1). 

Other kinds of inverses for difference operators are possible (see, e.g., [19J). The inverses 
are normally used in deriving hierarchies of soliton equations, in particular, non-isospectral 
hierarchies. 

Let u = u(n,t) = (ui(n, t), ■ • ■ ,u q (n,t)) T be a vector potential, in which n e Z and 
iGl are the lattice variable and the time variable, respectively. When an object P (e.g., a 
function or an operator) depends on u, its Gateaux derivative with respect to u in a direction 
v — (vi, ■ ■ • , v q ) T is defined by 

P'[v] = P'(u) [v] = —P(u + ev)\ £=Q = —P{ Ul + £Vl ,..-,u q + ev q )\ £=Q . (1.3) 

We denote by B the space of functions which are C°°-differentiable with respect to n and t 
and C°°-Gateaux differentiable with respect to u, and define the Lie bracket [•, •] on B q = 
{(Pi, ■ ■ • , P q ) T \Pi eB, 1 < i < q] as follows: 

[K, S] = K'[S] - S'[K] = K'{u)[S] - S'{u)[K], K,SeB q . (1.4) 

The forward difference operator D = E — 1 yields an equivalence relation ~ on B: 

P~Qif 3 ReB such that P-Q = DR. 

Let Xlnez P denote the equivalence class to which P belongs: 

J2 P = {P + DR | R e B}, PeB, (1.5) 

and J 7 , the quotient space: T = J-{B) = {J2 neZ P\P G £>}. An equivalence class of B by ~ 
is called a functional. The variational derivative 4? e B q of a functional V G JF with respect 
to m is determined by 



£(£)ViW)U^ 



It is easy to see that 



where P = P(u) e B. 
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^E^E^S' 1 ^^' (16) 
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The adjoint operator JT : B q —>■ B q of a linear operator J : B q — ► B q is determined by 

If Jt = —J, then J is called to be skew-symmetric. A linear skew-symmetric operator 
J : B q — > B q is called to be Hamiltonian, if the corresponding Poisson bracket 

{v , Q] = {v , eh = j : Q T j^ V]ef ,r, (i.7) 

satisfies the Jacobi identity: 

{V, {Q.,n}} + cycled, Q,K) = 0, P, Q,ft G T. 

A system of evolution equations u t = K, K G B q , is called to be a Hamiltonian system, if 
there are a Hamiltonian operator J : B q ^ B q and a functional 7i G JF, such that 

«t = AT = J-^. (1.8) 

The functional 7i is called a Hamiltonian functional of the system, and we say that the 
system possesses a Hamiltonian structure. 

We now assume that a pair of matrix discrete spectral problems 

E(f>=U(j) = U(u,\)(j), 

(1-9) 
(j) t = Vcj) = V(u, Eu, E- X u, • • • ; \)<f>, 

where u = u(n,t) is the potential, <p t denotes the derivative with respect to t, U, V G G are 
called a Lax pair and A is a spectral parameter, determines a discrete soliton equation 

u t = K = K(n, t, u, Eu, E- l u, ■■■), K eB q , (1.10) 

through their isospectral (i.e., A* = 0) compatibility condition (i.e., discrete zero curvature 
equation): 

U t = (EV)U-VU. (1.11) 

This means that a triple (U, V, K) satisfies 

U'[K] ={EV)U-VU, 

where £/"'[!£] denotes the Gateaux derivative of U with respect to u in a direction K. The Lie 
algebraic structure for such triples was discussed [18] and applied to non-isospectral flows 

EDI. 



To generate integrable couplings of the equation ( ll.lOp . take a semi-direct sum of G with 
another matrix loop algebra G c as introduced in [1]: 

G = G(tG c . (1.12) 



The notion of semi-direct sums implies that G and G c satisfy 

[G, G c \ C G c , 

where [G, G c ] = {[A, B] \ A E G, Be G c }. Obviously, G c is an ideal Lie sub-algebra of G. 
The subscript c indicates a contribution to the construction of integrable couplings. We also 
require that the closure property between G and G c under the matrix multiplication 

GG C , G C G C G c , 

where GiG 2 = {AB \ A G Gi, Be G 2 }, to guarantee that the discrete zero curvature 
equation over semi-direct sums of Lie algebras can generate coupling systems. 

Then choose a pair of enlarged matrix discrete spectral problems 

_ __ _ _ (1.13) 

(j) t = V(j) = V(u, Eu, E~ l u, • • ■ ; \)<f>, 

where the enlarged Lax pair is given by 

U = U + U C , V = V + V C , U C ,V C EG C . (1.14) 

Obviously, under the soliton equation (ll.lOp . the corresponding enlarged discrete zero cur- 
vature equation 

U t = (EV)U-VU (1.15) 

is equivalent to 

U t = (EV)U-VU, 

(1.16) 
U c>t = [(EV)U C - U C V] + [{EV C )U - UV C ] + \{EV C )U C - U C V C ). 

The first equation above precisely gives the equation fll.lOp . and thus, the whole system 
provides a coupling system for the equation f ll.lOp . This shows the procedure of generating 
discrete integrable couplings through semi-direct sums of Lie algebras, proposed in [3] . 

As usual, a bilinear form (-, ■) on a vector space is said to be non-degenerate when if 
(A, B) = for all vectors A, then B = 0, and if (A, B) = for all vectors B, then A = 0. Since 
semi-direct sums of Lie algebras are not semisimple, the Killing form is always degenerate 
on semi-direct sums of Lie algebras [12] , and thus, it is not helpful in analyzing Hamiltonian 
equations by the trace identity [2U E2] • Indeed, semi-direct sums of Lie algebras can carry 
particular algebraic structures [23], and the corresponding groups can extend the Poincare 
group to unite geometrical with internal symmetries in a nontrivial way [24]. A natural 
question here for us is whether we can replace the Killing forms with general bilinear forms 
to establish Hamiltonian structures for discrete soliton equations associated with semi-direct 
sums of Lie algebras. 

In this paper, we would like to answer this question. As in the case of the continu- 
ous variational identity [25], we are going to show that a discrete variational identity also 
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ubiquitously exists in discrete spectral problems and plays important roles in constructing 
Hamiltonian structures and thereby conserved quantities for discrete soliton equations. The 
crucial step of our success is that while presenting a discrete variational identity under a 
general bilinear form (•, •) on a given algebra g, we get rid of the invariance property 

(p(A),p(B)) = (A,B) (1.17) 

under an isomorphism p of the algebra g, but keep the symmetric property 

(A,B) = (B,A) (1.18) 

and the invariance property under the multiplication 

(A,BC) = (AB,C), (1.19) 

where AB denotes the product of A and B in g. If g is also associative, then g forms a Lie 
algebra under 

[A, B}= AB- BA, 

and the invariance property under the Lie bracket holds: 

(A,[B,C]) = ([A,B],C). (1.20) 

Conversely, the invariance property under the Lie bracket, ( 11.201) . doesn't imply the invari- 
ance property under the multiplication, ( 11.191) . We will show by examples that generally, 
there are many non-degenerate bilinear forms with the properties (11.181) and (11.191) on a 
given semi-direct sum of Lie algebras. 

The paper is organized as follows. First, in Section 2, we would like to establish a discrete 
variational identity under general non-degenerate, symmetric and invariant bilinear forms, 
in order to construct Hamiltonian structures of soliton equations associated with semi-direct 
sums of Lie algebras. Moreover, in Section 3, the constant 7 appeared in the variational 
identity will be determined precisely. Then, in Section 4, an application is given to a kind of 
semi-direct sums of Lie algebras in the Volterra lattice case, and consequently, Hamiltonian 
structures of the associated integrable couplings of the Volterra lattice hierarchy are pre- 
sented. This also justifies that the approach of integrable couplings using semi-direct sums 
of Lie algebras [I] can engender integrable Hamiltonian equations. A few concluding remarks 
are given in the last section. 

2 A discrete variational identity under non-Killing forms 

For a given spectral matrix U = U(u, A) e G, where G is a matrix loop algebra, let us fix 
the proper ranks rank(A) and rank(ii) so that U is homogeneous in rank, i.e., we can define 
rank(£7). The rank function satisfies 

rank(AB) = rank(A) + rank(S), 



whenever an expression AB makes sense, e.g., EU . Therefore, to keep the rank balance in 
equations, we have to define 

rank(E) = rank(C7) = 0. (2.1) 

The requirement rank(E') = is due to the stationary discrete zero curvature equation 

(EV)(EU) = UV, (2.2) 

and then the requirement rank(U) = is due to the discrete spectral problem E<fi = U4> in 

(H2||. 

Let us next assume that if two solutions V\ and Vi to ( 12.21) possess the same rank, then 
they are linearly dependent of each other: 

V 1 = 7 1/ 2 , 7 = const. (2.3) 

This is a strict condition on spectral problems, also required in deducing the trace iden- 
tity [21], the so-called quadratic- form identity [26] and the continuous variational identity 
[23] . which can be used to construct Hamiltonian structures of various continuous soliton 
equations (see, e.g, |27j-[29|). 

Associated with a non-degenerate bilinear form (-, •) on G with the symmetric property 
fll.ISp and the invariance property under the multiplication, (j!.19p . we introduce a functional 

W = X)«V; U x ) + (A, (EV)(EU) - UV)), (2.4) 

while U\ denotes the partial derivative with respect to A, and V,A £ G are two specific 
matrices. The variational derivative VaE G G of a functional TZ with respect to A e G is 
defined by 

Y(V A K,B) = -?-K(A + eB) , B e G. (2.5) 

*-^ OS £=0 

nez 

Obviously, based on the non-degenerate property of the bilinear form (•,-), we can have 

V B X)<A B) = A V B J^iA EB) = E- X A. 

It then follows from the symmetric property (11. IS)) and the invariance property under the 
multiplication, (11.191) . that 

VyW = U X + U(E- X A) - AU, V A W = (EV)(EU) - UV. (2.6) 

Note that the first variational derivative formula can not be obtained, if we only have the 
invariance property under the Lie bracket, ( 11.201) . 



2.1 A discrete variational identity 

We are going to prove that there is a variational identity in the discrete world, similar to 
the continuous variational identity [25] . 

Theorem 2.1 (The discrete variational identity under general bilinear forms): Let G be a 
matrix loop algebra, and U = U(u, A) G G be homogenous in rank such that ^2.2\) has a 
unique solution V G G of a fixed rank up to a constant multiplier. Then for any solution 
V G G of $2.2\) . being homogenous in rank, and any non- degenerate bilinear form (•, •) on 
G with the symmetric property (11.181) and the invariance property under the multiplication, 
( I1.19p . we have the following discrete variational identity 

where -f- is the variational derivative with respect to the potential u and 7 is a constant. 

Proof. Let us start with the functional W introduced in ( 12. 4ft . For the variational calculation 
of W with respect to the potential u, we require the following constraint conditions: 

V V W = U x + U(E-^k) - AU = 0, (2.8) 

V A W=(EV)(EU)-UV = 0, (2.9) 

to determine V and A. These conditions also imply that both V and A are related to U and 
thus to the potential u. Immediately from the second constraint condition (I2.9p . we have 

Now using both of the constraint conditions (12.81) and (12.91) . and noting the property that 
if V aR.(A) = 0, then j-1Z(A(u)) = for a functional 7£, we know that only the dependence 
of u in U (but not in V and A) needs to be considered in computing -^-. Therefore, based on 
the invariance property under the multiplication, (I1.19P (note that the invariance property 
under the Lie bracket, ( ll.20p . is not good enough), we obtain 

'^^ » = ( ^ ) + ( e,5, (2.10) 

ou ^— ' ou ou ou 

where 

6 = (E- l A)V - VA. (2.11) 

This matrix satisfies 

E(e-v x )(EU)-u(e-v x ) = o, 



namely, Q — V\ solves (12.21) . This is because we have 

(EQ)(EU) - U& = A(EV)(EU) - (EV)(EA)(EU) - U{E~ 1 A)V + UVK 

= AUV - (EV)(EA)(EU) - U{E- X A)V + (EV)(EU)A 

= [AU - UiE^A^V + (EV)[(EU)A - (EA)(EU)} 

= U X V - (EV)(EU X ) 

from (J23D and $ZM), and 

(EV X )(EU) - UV X = U X V - (EV)(EU X ) 

from differentiating ( 12.91) with respect to A. By taking use of the uniqueness condition (12.31) 
and rank(0 — VA)=rank(VA)=rank(^l / ), there exists a constant 7 such that 

Q-V x = (E~ 1 A)V-VA-V X = ?-V, (2.12) 

A 

because jV is also a solution to (12. 2ft . 

Finally, ( 12.101) can further be expressed as 

This completes the proof. □ 



2.2 A formula for the constant 7 

Let us consider the other form of the stationary discrete zero curvature equation 

(ET)U-UT = 0. (2.13) 

If V is a solution to (12. 2p . then T = VU satisfies 

DT = [U,V], (2.14) 

where D = E — 1, as defined in the introduction. This is a counterpart of the stationary 
continuous zero curvature equation V x = [U, V}. When U is invertible, then V is a solution 
to (D iff T = VU is a solution to (l2~T3l . 

The matrix V^ presents the gradient which is needed for constructing the desired Hamil- 
tonian structure, and the matrix V contributes to the constant 7 in the variational identity 
as follows. 



Theorem 2.2 Let V be a solution to ( \2.2\i and Y = VU. Then for any bilinear form (-,-) 
on G with the properties ( II. 18ft and ( 11. 20ft . we have 

D(T m ,T m ) = (E-l)(Y m ,Y m ) = 0, m>\. (2.15) 

Proof. Noting that EY = E(VU) = UV, it follows from the symmetric property (11.181) and 
the invariance property under Lie bracket, ( 11.201) . that 

D(Y m ,Y m ) = ({UV) m ,{UV) m )-({VU) m ,{VU) m } 

= {{UV) m -{VU) m ,{UV) m + {VU) m ) 

= ([U,V{UV) m -\{UV) m + {VU) m ) 

= {U, [v{uv) m -\ {UV) m + {VU) m }) 

= (U, [v{uv) m - x {vu) m - x v, U\) 

= (ViUV^iVU^V, [U, U]) = 0, 
where m > 1. This proves the theorem. □ 

By ( 12.151) . (r, T) is independent of the lattice variable n. 

Theorem 2.3 Let V be a solution to (12.21) and Y = VU. If (Y, Y) ^ 0, then the constant 7 
in the discrete variational identity ( 12.71) is given by 

7=-^ln|<r,r>|. (2.16) 

Proof. It follows from (EH]) and (l2~T2l) that 

r A = (vu) x = v x u + vu x 

= [(^ 1 A)F - VA - ^-V}U + V[MJ - U(E- X K)} 

A 

= [E- 1 \,VU]-'\VU 
A 

= [E- 1 A,Y]-^Y. 
Therefore, differentiating (r,T) with respect to A yields 

(r,r> A = (r A ,r> + (r,r A > = 2(r A ,r> 

^2([E- 1 A,Y]-^Y,Y) 

= 2([E- 1 A,Y],Y)-^(Y,Y) 
= 2(E- 1 A,[Y,Y])-^(Y,Y) 

= -|<r,r>. 

This implies that (12~TB|) holds. □ 

Note that the formula ( 12.161) for the constant 7 in (" 12. 121) is still true, if we only have the 
invariance property under the Lie bracket, (11.201) . 
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3 Symmetric and invariant bilinear forms 

Let us consider the following semi-direct sum of Lie algebras of 4 x 4 matrices 



G = G E G r 



A 
A 





A = 


as a 4 


}•{ 


A x 



A 1 = 


a 5 a 6 
a 7 a 8 



,(3.1) 



where Oj, 1 < i < 8, are real constants. In order to construct symmetric and invariant 
bilinear forms on G conveniently, we transform the semi-direct sum G into a vector form. 
Define the mapping 



a : G -> R s , A i-» (oi, • • • , a 8 ) T , A 



d\ (l 2 Q>5 <^6 

0-3 0,4 Qj-j tig 

a\ a 2 

03 a 4 



EG. 



(3.2) 



This mapping a induces a Lie algebraic structure on R 8 , isomorphic to the matrix loop 
algebra G. The corresponding Lie bracket [■, ■] on M 8 can be computed as follows 



[a, b] T = a T R(b), a = (oi, • • • , a 8 ) T , 6 = (b u ■ ■ ■ , b 8 ) T E M 8 , 



(3.3) 



where 



R(b) 



6 2 -63 b 6 -b 7 

63 b^ — b\ —63 67 b$ — b 5 —67 

-62 61 — 64 b 2 —b 6 h- b$ b e 

-62 63 -6 6 67 

b 2 -63 

63 64 - h -63 

-b 2 61 - 64 62 

-62 63 

This Lie algebra (1R 8 , [•,•]) is isomorphic to the matrix Lie algebra G, and the mapping a, 
defined by f 1 3 . 2 j) . is a Lie isomorphism between the two Lie algebras. 

A bilinear form on M. 8 can be defined by 

(a, b) = a T Fb, (3.4) 

where F is a constant matrix (actually, F = ((ej, e,,)) 8x8 , where e 1 , ■ ■ ■ , e 8 are the standard 
basis of M 8 ). The symmetric property (a, b) = (b, a) requires that 



F T = F. 



(3.5) 



Under this symmetric condition, the invariance property under the Lie bracket 

(a, [b,c]) = ([a,b],c) 
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equivalently requires that 

F{R{b)) T = -R{b)F, b e R 8 . (3.6) 

This matrix equation leads to a linear system of equations on the elements of F. Solving the 
resulting system yields 



m o 



r/i — 772 



77 2 Vi 

771-772 



V'2 



















771 774 

V<i V5 



773-774 



Vi 







773-774 


773 775 














774 

773-774 













V5 




V5 



(3.7) 



where 77^, 1 < i < 5, are arbitrary constants. Now, the corresponding bilinear form on the 
semi-direct sum G of Lie algebras defined by (13.1 p is given as follows 



(A,B) 



G 



(a-\A),a- 1 (B)) R s = (a ir ■ ■ ,a 8 )F(b u - ■ ■ ,b 8 ) T 

(77^1 + 77 2 a 4 + 77305 + 77 4 a 8 ) 61 + [(771 - 772) a 3 + (773 - 774) a 7 ] b 2 

+ [{Vi - V2) «2 + (V3 ~ Vi) Oe] h + (77204 + 771 a 4 + 7? 4 a 5 + ^s) h 

+ (%ai + Vi a 4 + %«5 + V^as) h + (r/3 ~ Vi) «3^6 

+ (V3 ~ Vi) a 2 b 7 + (r?4Gi + 773^4 + %a 5 + %«s) b 8 , 



(3. 



where 



A 



fl-l 


«2 


a 5 


a 6 




&i 


lh 


&5 k 


a 3 


°i 


a- 


a 8 


, B = 


h 


h 


b 7 b 8 








Oi 


a 2 










bi b 2 








03 


04 










b 3 b 4 



EG. 



The bilinear form (13.81) is symmetric and invariant under the Lie bracket of the matrix 
Lie algebra: 

(A,B) = (B,A), (A,[B,C\) = ([A,B],C), A,B,CeG. 

But this kind of bilinear forms is not of Killing type, since the matrix Lie algebra G is not 
semisimple. A direct computation shows that the bilinear form (I3.8P is also invariant under 
the matrix multiplication: 

(A, BC) = (AB, C), A,B,Ce G. 

We started with the invariance property under the Lie bracket but not under the multiplica- 
tion, since it is easier to express the invariance property under the Lie bracket as an equation 
like fl3~6|) . 
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The bilinear forms defined by (13.81) contain plenty of non-degenerate cases. A particular 
non-degenerate bilinear form with rji = 772 = ^3 = 1 and 774 = 775 = will be used to 
establish Hamiltonian structures for the integrable couplings of the Volterra lattice hierarchy 
associated with the above semi-direct sum of Lie algebras. 



4 Application to the Volterra lattice hierarchy 



4.1 The Volterra lattice hierarchy 



Let us recall the Volterra lattice hierarchy 
Volterra lattice hierarchy is given by 

E<j)=U(f), U = U(u,X) 



A discrete spectral problem for the 



1 u 
A- 1 



This is equivalent to 
Upon setting 



\{E 2 - E)(f) 2 = uh- 



a b 
c —a 



E^-=E 



«>o 



j>0 



Ca (la 



\- 



(4.1) 



(4.2) 



the discrete stationary zero curvature equation (12.131) gives rise to 

' a® + A -1 &W - a - uc = 0, 

Ma' 1 ' — b + ua = 0, 
C (D _ A-V 1 ) - \-'a = 0, 
^ mc« - X^b = 0, 
which equivalently leads to 

b = u(a {1) + a), c = A~ x (a + o (_1) ), 

a« - a + \- l [u^(a^ + a«) - u(a + a^)] = 0. 

This system can uniquely determine all sets of functions Oj, b{ and q, upon choosing 

1 
oo = ■x, ai\u=o = 0, 2 > 1. 

In particular, the first two sets are 

oo = |, b = u, c = 0; 



(4.3) 



(4.4) 



a 1 



-u, 61 



u(u^ + u), C\ = 1. 
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The compatibility conditions of the matrix discrete spectral problems 



E<p = U<f>, & = ^ [m V, V^ = (A m+1 r) + + A m , A r 



—bm+1 

O-m+1 + a m+l 



m > 0, (4.5) 



where (A m+1 r) + denotes the polynomial part of A m+1 T in A, determine (see, e.g., [18]) the 
Volterra lattice hierarchy of soliton equations 



u tm =K m = ^ m K = u(a% +1 - a { -l\), K = u(u^ - u«), m > 0, 
where the hereditary recursion operator $ is given by 

$ = w(l + E- l ){-u {x) E 2 + u){E - l) -1 !* -1 . 
Since we have 

(V, U x ) = tr(VU x ) = \- l a^, (V, U u ) = ti{VU u ) = ~ 



(4.6) 



(4.7) 



u 



where V = TU' 1 , an application of the trace identity with 7 = in [21] (corresponding to a 
particular case of (13.81) : 771 = 1 and rji — 0, 2 < i < 5) presents the Hamiltonian structures 
for the Volterra lattice hierarchy: 



sn 



u tm = J^, J = u(E-> - E)u, H m = D-^tH m * °" 



5m 



n€ < 



m + 1 



(4.8) 



4.2 Hierarchy of integrable couplings and its Hamiltonian struc- 
ture 



As in pEJ, introduce two Lie algebras of 4 x 4 matrices 

AeR[A]®gl(2)l, G c = 



G 



A 
A 



B 




B e R[X] ® gl(2) 



(4.9) 



where the loop algebra R[A] ® gl(2) is defined by span{A n /L \n > 0, A e gl(2)}, and form 
a semi-direct sum G = G (E G c of these two Lie algebras G and G c . In this case, G c is 
an Abelian ideal of G. For the Volterra spectral problem (14.11) . we define the corresponding 
enlarged spectral matrix as follows 



U=U{u,\) 



u u a 

u 



EG(±G C , U a = UJv) 



v 




(4.10) 



where v is a new dependent variable and the enlarged potential u reads 

u = (u, v) T . 



(4.11) 
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To solve the corresponding enlarged stationary discrete zero curvature equation 

(ET)U-tJf = 0, (4.12) 

we set 



r r 
o r 



, r a = r a (w, A) 



e / 



9 ~e 



(4.13) 



where T is a solution to (12.131) . defined by (14. 2p . Then, the enlarged stationary discrete zero 
curvature equation gives 

[(ET a )U - UT a ] + [(ET)U a - U a T] = 0, 

together with (12.131) . This equation equivalently leads to 

' eW + A -1 / (1) - e - ug - vc = 0, 
v aP~' + ue^ — f + ue + va = 0, 
g (i) _ A -i e (i) _ A -i e = o, 
vcW + ug 0) - A" 1 / = 0. 

Since the second equation is always satisfied if the last two equations hold, this system is 
consistent and determines 



Trying a solution 



and choosing 



/ = u ( e (i) + e ) + v(aW + a), g = \-\e + e^ 1 )), 

e (D _ e + A- 1 [ M «(e( 2 ) + e«) + ^«(a( 2 ) + aP>)\ - uc = 0. 

e = J2 e * A " f = E ^ A " » = E ^ A "> 



i>0 



i>0 



j>0 



e = 0, ei\ Q=0 = 0, i > 1, 



(4.14) 



(4.15) 



(4.16) 



we see that all sets of functions e^, ft and gi are uniquely determined. In particular, the first 
two sets are 

e = 0, /o = v, g = 0; 

ei = — 1>, /i = — m(v^^ + v ) — t>('u( 1 ) + u), g\ = 0. 



Let us now define 



V'l' 



W n 



e G, K H = (A m+1 r Q ) + + A m , a , m > 0, 



(4.17) 



where y' m ' is defined as in (14.51) and (\ m+1 T a ) + denotes the polynomial part of X m+l Y a in 
A, and choose A m>a as 

—fm+i 



<-^m,a 



U e m +i + e m +l 



m > 0. 



(4.18) 
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Then, the m-th enlarged discrete zero curvature equation 

U tm = (EV [m] )U - UV [m 
leads to 



vu 



,(-i) 



(-1). 



/ m +i - u(e m+1 + e m+1 ) - v(a m+1 + a m+1 ) 



,W 



,(-!) 



(1) 



,(-!)• 



M ( e m+1 e m+l) + u ( a m+l a m+l) 



(4.19) 



together with the m-th Volterra lattice equation in (14. 61) . Here in the last equality, we used 
f)4.14p . A hierarchy of coupling systems are thus generated for the Volterra lattice hierarchy 



Ut r , 



~ 




= K m (u) = 


= $ m K = 







(-i)n 

M l a m+1 a m+l/ 



u(e. 



(i) 



"m-t 



m+1 C m+L 



,W 



(-I)- 



+ U ( a m+1 a m+l, 



171 > 0, 



(4.20) 



in which the first system w io = i^o reads 



M<o = u ( u (-i) - u (i)), Uto = v ( u (-i) - u (i)) + „(„(-!) - „(i); 



and the hereditary recursion operator $ is defined by 

$ 



$ 



$„ $ 



(4.21; 



where $ is given by ( 14. 7ft and 

$ c = u(l + r')(-t( (1) £ 2 + u)(E- l)- 1 ^ 1 + w(l + E'^i-v^E 2 + v)(E- l)" 1 -^ 1 
-w(l + E- l )(-u^E 2 + w)(£ - l)- 1 ^- 2 . 

To construct Hamiltonian structures of these integrable couplings by using the discrete 
variational identity (12.7J) . we consider a non-degenerate bilinear form on G = G € G c defined 
by (13.81) under the selection of 



m = V2 = V3 = 1 3 Vi = V5 = o. 
Then, a direct computation tells 



(V,Ux) 
(V,U U ) 



vb + u 2 e — uf 



A-V 1 ), 



;/- 



Am 2 
--, (F,C7„) =■■---. 



u 



w 



(4.22) 

(4.23) 
(4.24) 



where U is defined by (I4.10|) and V = TU x with F being defined by (14.131) . Now an appli- 
cation of the discrete variational identity ( 12.7ft with 7 = engenders 



8u 



Tir. 



'VClm+i e m _|_i 0"m+l\T 



u- 



tl 



U 



1 "ra 



Ei 



n€ i 



67714- 1 

m + 1 



j, m> 0. 



(4.25) 
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Consequently, we obtain the Hamiltonian structure for the hierarchy of integrable couplings 
in (lOnjl : 

(4.26) 



'■■■« 



ou 



where the Hamiltonian functionals TC m , m > 0, are given in (14. 25ft and the Hamiltonian 
operator J is determined by 



J 



^e^-e^ 

utE- 1 - E)u u(£ _1 - E)v + v{E~ l - E)u 



(4.27) 



Now, noting that $J = J$t, it follows that each Hamiltonian coupling system in the above 
hierarchy (14.261) possesses infinitely many conserved functionals {'H n }'^ =Q and infinitely many 
symmetries {K n }™ =0 , which commute with each other: 



{H k ,Ht} = 0, [#*,#,] = 0, k,l>0. 



(4.28) 



5 Concluding remarks 



The trace identity has been generalized to discrete spectral problems associated with non- 
semisimple Lie algebras, or equivalently, Lie algebras possessing degenerate Killing forms. 
The constant 7 in the discrete variational identity was determined precisely by the corre- 
sponding solution to the stationary discrete zero curvature equation. The resulting discrete 
variational identity was applied to a class of semi-direct sums of Lie algebras in the Volterra 
lattice case and furnished Hamiltonian structures for the associated integrable couplings of 
the Volterra lattice hierarchy. 

The proof of the variational identity involves, in an essential way, the invariance property 
of bilinear forms. We remark that there is a difference between the continuous and discrete 
cases. In the continuous case, we only require the invariance property of bilinear forms under 
the Lie bracket of the underlying algebras, but in the discrete case, we require the invariance 
property of bilinear forms under the multiplication. In theory, the invariance property under 
the multiplication is stronger than the invariance property under the Lie bracket, because 



(A,[B 1 C]) = ([A 1 B],C) 



is just a consequence of 



(A,BC) = (AB,C). 



One such example is the Lie algebra 



G 



Oi a 2 
a 3 



Qi G 



1< i < 3 
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for which the corresponding matrix F is 


V2 773 

where r/j, 1 < z < 3, arbitrary constants. Though all associated bilinear forms on this three 
dimensional Lie algebra 



(A, B) = (77101 + 772^3)61 + (V20-2 + V3 a s)h, A 



are invariant under the Lie bracket of matrices, they are varied under the matrix multiplica- 
tion. Nevertheless, we see that both invariance properties are equivalent to each other in the 
case of the semi-direct sum of matrix Lie algebras in Section [SI We also remark that more 
chooses of non-degenerate bilinear forms in ( 13.81) could lead to more Hamiltonian integrable 
couplings for the Volterra lattice hierarchy, and similar applications could be made for other 
lattice hierarchies 



Oi 


02 


, B = 


b\ b 2 





% 




h 



To conclude, the discrete variational identity ubiquitously exists in discrete spectral 
problems associated with both semisimple Lie algebras and non-semisimple Lie algebras. It 
brings us a powerful tool for exploring Hamiltonian structures of discrete soliton equations. 
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